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RESUMEN 

En este trabajo evaluamos la masa mmima de ciimulos estelares, ^Vf™™^ para la cual los efectos de muestreo de 
la Funcion Inicial de Masa (FIM) no pueden ignorarse. Esta tVI™*" corresponde a la situacion para la cual la 
luminosidad del cumulo es igual a la luminosidad de la estrella individual mas brillante incluida en el modelo, 
y toma valores entre 3 x 10^ y 6 x 10^ M© segiin la edad y la banda observada. Mostramos ejemplos para 
poblacioncs jovenes (i < 10 Ma) y viejas {t hasta 10 Ga). Tambien hacemos publica la distribucion espectral 
de energia (DEE) desde 160 /^m (1.8 x 10^^ Hz) hasta 25 keV (6.0 x 10^* Hz) para regiones de formacion estelar 
con edades entre 0.1 y 10 Ma y mctalicidades entre Z=0.001 y ^=0.040. Las DEEs estan accesibles en nuestro 
servidor WWW e incluyen las cantidades necesarias para cvaluar los efectos de muestreo de la FIM. Las DEEs 
pueden usarse tanto para la obtencion de colores (incluyendo los de rayos X), como para entrada de modelos 
de fotoionizacion o para ajustes de datos observados tcnicndo en cuenta la incertidumbrc intrmseca en los 
modelos debida al muestreo de la FIM. 

ABSTRACT 

In this contribution we evaluate the minimal cluster mass for which the effects of an incomplete sampling of the 
Initial Mass Function (IMF) cannot be neglected. This minimal cluster mass corresponds to the situation in 
which the integrated luminosity of the cluster modeled equals the luminosity of the most luminous individual 
star included in the model, and it takes values between 3 x 10^ and 6 x 10^ M0 depending on the age and the 
observed band. We show different examples for young {t < 10 Myr) and old {t up to 10 Gyr) stellar populations. 
We also make here a first release of the spectral energy distribution (SED) from 160 (1.8 x 10^^ Hz) to 
25 keV (6.0 x 10^® Hz) for star-forming regions with ages between 0.1 and 10 Myr, and metallicities between 
Z=0.001 and Z=0.040. The SEDs are available at our WWW server and include the corresponding quantities 
for the evaluation of sampling effects. These SEDs can be directly used to obtain colors (including the X-rays 
ones), as input of photoionization codes, or for fitting with observed data taking into account the intrinsic 
uncertainty of the models due to the IMF sampling. 
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1. THE NEED FOR A STATISTICAL 
MODELING 

Evolutionary synthesis models have been exten- 
sively used since the work of Tinsley & Gunn (1976). 
Unfortunately, the limitations for their use have not 
been extensively studied. In this work we present 
the most basic limit for the use of standard synthesis 
models (those that assume a continuously populated 
Initial Mass Function, IMF). A more detailed study 
can be found in Cerviho & Luridiana 2003 (CL03, 
submitted) . 

Maybe, the most trivial limit for the usage of a 
synthesis model is the following one: The total lumi- 
nosity of the cluster modeled must be larger than the 
individual contribution of any of the stars included in 
the model. This obvious statement defines a natural 
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theoretical limit that has not always been considered 
when synthesis models are applied to real observa- 
tions. Based on this limitation, we can establish a 
Lowest Luminosity Limit (LLL) for the application 
of synthesis models, which corresponds to the situ- 
ation where the integrated luminosity of the cluster 
modeled equals the luminosity of the most luminous 
individual star included in the model. 

Since the integrated luminosity scales with the 
initial amount of gas transformed into stars, M, a 
minimal initial amount of gas transformed into stars, 
A^™™, can be inferred for the use of synthesis mod- 
els. In clusters with masses below Al™*", sampling 
effects in the IMF cannot be neglected and a statisti- 
cal modeling is needed. Even more, it can be demon- 
strated that even the results of synthesis models for 
clusters with < 10 x A^™*" have an intrinsic rel- 
ative uncertainty equal or larger than 10%. 
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Fig. 1. TVf""" as a function of the age for different pho- 
tometric bands. 

2. SOME EXAMPLES OF THE LLL 

In order to illustrate how relevant this efFcct may 
be, we show in Fig. ^the values at different 

ages and photometric bands. In this case the ^Vf™*" 
values have been computed from the isochrones by 
Girardi et al. (2002) with metaUicity Z=0.0004, and 
simple stellar populations (or Instantaneous Burst, 
IB) results, and it is assumed a Kroupa (2001) IMF 
in the mass range 0.01 - 120 Mq. 

In the case of young stellar populations (ages 
lower than 10 Myr) we show in Fig. [2 the X™™ 
values corresponding to the number of ionizing pho- 
tons above the H°, He" and He+ ionization edges. 
The tVC"™ values correspond to models following a 
Salpetcr (1995) IMF in the mass range 0.1 - 120 Mq 
and they have been computed from a hybrid version 
of the codes by Cerviho, Mas-Hesse & Kunth (2002, 
CMHK02) and Leitherer et al.(1999) in order to both 
compute yV^™™ and include the atmosphere models 
from Smith, Norris, & Crowther (2002). In the fig- 
ure we use evolutionary tracks with high mass-loss 
rates from Meynet et al. (1994). X-ray emission has 
not been considered in this figure. 

Note that, in the galactic context, the most 
massive OB association known, Cygnus 0B2 (c.f. 
Knodlseder 2000), has transformed into stars ~ 10^ 
Mq of gas. Hence, the modeling of normal galactic 
and extragalactic HIT regions, with a lower amount 
of gas transformed into stars, is potentially affected 
by sampling effects. 

3. STATISTICAL SYNTHESIS MODELS 

As we have seen before, the IMF sampling may 
play an important role in the determination of the 
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Fig. 2. Al™'" as a function of the age for ionizing photons 
above the H'', He'' and He+ edges. 

evolutionary status of stellar clusters obtained via 
synthesis models. Some studies from a qualitative 
point of view can be found in, e.g., Cervifio, Lurid- 
iana, & Castander (2000), Bruzual (2002), and ref- 
erences therein. However, a quantitative theoretical 
formalism is needed to address this subject and ap- 
ply the results to real data. 

Examples of such kind of formalism can be found 
in Langon & Mouhcine (2000), or in Cerviho et al. 
(2002). These last authors, based on the definition 
by Buzzoni (1989) of an effective number of stars, J\f, 
evaluate quantitatively the sampling effects in stel- 
lar clusters. The formalism is completely valid for 
any age range and for quantities that scale linearly 
with A4 . It also gives a first order estimation of the 
bias in the ratios or logarithmic quantities predicted 
by synthesis models when the sampling effects are 
quite important (Cerviho & Valls-Gabaud 2003), or, 
cquivalently, when the observed cluster have A4 val- 
ues close to the LLL (see CL03 for more details). 

Now, we make here a first public release of the 
complete spectral energy distribution (SED) and 
the corresponding Af values at each wavelength 
needed to evaluate the sampling effects. The re- 
sults come from an improved and extended version 
of the synthesis code presented in Cervifio & Mas- 
Hesse (1994), and have been discussed in CMHK02. 
A first experimental public release of the code will 
be available in the following months. 

The SEDs range from 160 /xm (1.8 x lO^^ Hz) 
to 25 keV (6.0 x 10^^ Hz), for ages between 0.1 and 
10 Myr and for metallicities between Z=0.001 and 
Z=0.040. They can be directly used to obtain colors 
(including the X-rays ones), as input of photoion- 
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Fig. 3. Multiwavelength emission for a 3.2 Myr star- 
forming region with the corresponding 90% Confidence 
Limits. A bremsstrahlung component with 10*^ K and 
a total luminosity of 4 x lO'*" erg s~^ is also shown for 
comparison. 

ization codes, or for fitting with observed data 
taking into account the intrinsic uncertainty of the 
models due to the IMF samphng. The data, together 
with the samphng effects for other observables, can 
be obtained in tabular form at: 

|http : //www. laef f .esa. es/users/mcs/SED/| 

One of the most important features of our newly 
realesed SEDs is the extension in the X-ray do- 
main. As an example of a possible application of 
this feature, we considered the work by Stasihska 
& Izotov (2003), who study the conditions needed 
to reproduce the observed sequences of low metal- 
licity star- forming regions, address the difficulties 
that traditional models (the ones without X-rays) 
would find to reproduce at the same time all the 
observational constraints, and point out directions 
where one should go in order to find physical solu- 
tions to the problems encountered. Among other ef- 
fects, they show that an additional high-energy com- 
ponent is needed to explain the observed trend in 
diagnostic diagrams, and, as an example, they use 
a bremsstrahlung spectrum with 10® K and a to- 
tal luminosity of 4 x 10'*° erg s~^ (Fig. |2| dashed 
line). Incidentally, this component corresponds ac- 
tually to the type of supersoft X-ray sources ob- 
served by ROSAT (Rappaport et al. 1994). In Fig. 
|31 we show the multiwavelength emission predicted 
by our code for a 3.2 Myr star-forming region with 
Z^O.OOl and X=2.3xl0^ Mq (i.e., M=10''' Mq in 
the mass range 0.8 - 120 Mq, as in Stasihska & Izo- 
tov 2003) for different efficiencies of conversion of ki- 
netic energy in X-rays. The kinetic energy included 



in our model is produced by stellar winds and su- 
pernovae (CMHK02), so it is totally self-consistent 
with the stellar population responsible for the ion- 
ization of the gas. The comparison between the 
bremsstrahlung component by Stasihska & Izotov 
(2003) and the high-energy portion of our spectrum 
shows that the availability of SEDs extending into 
the high-energy range might help to explain some 
observed features of stellar populations; more specif- 
ically, the newly released SEDs will give the opportu- 
nity to perform a more detailed analysis of the effects 
of X-rays in the optical emission line spectrum. 
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